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ABSTRACT 

Planetary migration poses a serious challenge to theories of planet formation. In 
gaseous and planetesimal disks, migration can remove planets as quickly as they form. 
To explore migration in a planetesimal disk, we combine analytic and numerical ap- 
proaches. After deriving general analytic migration rates for isolated planets, we use 
A^-body simulations to confirm these results for fast and slow migration modes. Mi- 
gration rates scale as m"^ (for massive planets) and (1 -|- (ch/B)^)"^, where m is the 
mass of a planet and en is the eccentricity of the background planetesimals in Hill units. 
When multiple planets stir the disk, our simulations yield the new result that large-scale 
migration ceases. Thus, growing planets do not migrate through planetesimal disks. To 
extend these results to migration in gaseous disks, we compare physical interactions 
and rates. Although migration through a gaseous disk is an important issue for the 
formation of gas giants, we conclude that migration has little impact on the formation 
of terrestrial planets. 

Subject headings: planetary systems - solar system: formation - stars: formation - 
circumstellar matter 



1. Introduction 



Migration is an important physical process in planet formation (e. 



Lin &: Papaloizou 



Kirsh et al 



1986; 



2nn9l : lDAngelo et al. 



Ward I1997I: lArtvmowiczjbnni iLevison et allbnnTi : IPapaloizou et al.lbnOTl : 

20ld : iLubow k Idalboiol . and references therein). Based on analytic theory and detailed numerical 
simulations, several modes of interaction between growing planets and density perturbations within 
a disk of gas or within a disk of planetesimals produce secular evolution of the orbital semimajor 
axis, eccentricity, and inclination of a planet. For planets with masses exceeding ~ 0.1 M0, derived 
migration rates have a broad range, ~ 10"'' — 10~^ AU yr~^. On typical timescales of 0.1-1 Myr, 
planets can migrate through the entire disk. 
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To explain the frequency of ice giant and gas giant planets c l ose to their parent stars, mi- 



grati on is essential (|Lin et alJll996l : lMarzari &: Weidenschillingj|2002l : llda &: Linll2004l : lAlibert et al 
2004 ). Although there are significant selection biases, most known exoplanets have semimajor 
axes, a < 0.1-1 AU (data from exopIanet.org and exoplanet.eu). Protostellar disks probably do 
not haye enough rnass t o produce ice gian ts or gas giants so close to their parent stars (e.g., 
Bodenheimer et al.ll200Cl : iKornet et al.ll2002l ) . Once these planets form farther out in the disk, how- 
ever, they can slowly migrate inward to close-in orbits around th eir parent stars (e.g.. llda fc: Lin 



2005 



Armitagd 120071 : iThommes et al 



2008 



Mordasini et al.ll2009l ). 



Migration may also explain the orbital architecture of the Solar System. Observations of the 
dynamical structure of the Kuiper belt suggest tha t Neptune migrated outward from its likely 
birthplace (IMalhotra I Il993l : iHahn &: Malhotral Il999l ) . Other evidence suggests t hat the four gas 



giants formed in a inore compact configuration a nd then migrated outward (e.g., iThommes et al 



2002 



Tsiganis et al. 



2005 



Morbidelh et al.ll2008l ). 



Despite these successes, migration is a great challenge for theories of planet formation. In 
the current picture, terrestrial planets and the cores of at least some gas giant planets form by 
a coagulation process, where lower mass objects collide and merge into larger objects. Early on, 
migration timescales are long. Without straying too far from their birthplaces, protoplanets undergo 
runaway growth - where a few of the largest objects grow much much faster than other objects - and 
then oligarchic growt h - where these larges t objects grow more slowly but still f aster than much less 



mass ive objects (e.g., iKokubo Ida 



1998 



Goldreich. Lithwick. &: Sari 



2004 



Kenyon Bromley 



2010l . and references therein). As planets begin to reach masses of ~ 0.1 M^, however, collision times 
become longer than migration times. Thus, theory predicts that the final building blocks of planets 



migrate into the central star before they reach the mass of the Earth (ILin &: Papaloizou 



197S; 



Goldreich fc Tremainelll98d : lArtymowic2]ll993bl : IWardlll997l : iMasset Papaloizoull2003l : llda Lin 
200 



Migration is also a severe problem for the formation of ice giant and gas giant planets. Once 
ice giants or gas giant s are fully- forme d, migration can produce the close-in giant planets observed 
around nearby stars (jlda Linll2005l ). Howev er, theory predicts a more rapid migration of the 



lower mass building blocks o f ice and gas giants (jWard 



2008 



1997 



Masset &: Papaloizoull2003l : llda &: Lin 



Peplihski et al.ll2008al ). In the standard theory, these lower mass planets migrate too fast to 



produce ice or gas giants. Solving this problem is a central issue in planet formation theories. 



Theories of migration generally focus on isolated planets interacting with the disk (see lPapaloizou et al 



20071 . and references therein) . Recent attempts to understand how real planets avoid migratio n con- 



centrate on the physics of this isolated interaction, inclu ding disk dynamics (iMasset et al 



Peplihski et al. 



2008a 



Paardekooper &: Papaloizou 



orbital eccentricity (iPapaloizou &: Larwoodll200o[K disk therinodynamics (iKley Sz Cridall^uugl. 
Paardekooper &: Mellema|l2006bl .|Paarde kooper fc Papaloizoul[20c3. IkIcv et ahlbood. Paardekoop 



2009a 



2009b), magnet ic fields (iTerquem 



2006,: 



200. 



2008, 



3), 



2010 . Paardekooper et al.boil ). and turbulence ( Nelson k Papaloizoj2004 : Adams Bloch 2009). 



)er et al 



While any or all of these processes may reduce migration rates to acceptable levels, growing pro- 
toplanets are not isolated. Tightly packed protoplanets probably perturb the disk differently than 
systems of widely spaced protoplanets. Thus, migration rates may depend as much on the local 



density of protoplanets as on the scale of specific interactions between an isolated planet and the 
disk. 

Here, we co nsider how migr a tion operates in systern s of multiple planets. Building on pre - 
vious work (e.g., Malhotra Ill993l : iHahn k MalhotralllQod : iLevison et ahlboOTl : Kirsh et ahlbood ). 
we examine migration in disks of planetesimals with single planets (§2.1-2.3) and multiple plan- 
ets (§2.4-§2.5). These results show that migration is rarely important in planet-forming disks of 
planetesimals. In §3, we then explore the implications of our results for (inviscid) planetesimals 
embedded in (viscous) gaseous disks. If our assumptions about viscous disks are valid, migration is 
rarely important during terrestrial planet formation. However, it is still an important issue in the 
formation of ice giant and gas giant planets. We conclude with a brief summary and suggestions 
for further study in §4. 



2. Planetary migration in a planetesimal disk 

Planets migrate through a planetesimal disk as a result of pairwise exchange of angular momen- 



1979, 


1980: 


Artvmowic2 


1993a^ 



radius. 



\3MJ ' 

where a is the planet's semimajor axis, m is its mass, and M is the mass of the central star. If the 
semimajor axis of a disk particle is a + dr, where 6r is its orbital separation from the planet, then 
a passing encounter changes the planet's semimajor axis by 



6d 



9[.x) — , 
m 



(2) 



where x = Sr/rn, g{x) is a function that depends on the geometric shape of the planetesimal's 
trajectory relative to the planet, and the tilde symbol indicates a change in orbital distance per 
unit planetesimal mass. Equation ([2]) asserts that the dimensions of a particle's trajectory near the 
planet scale as rn; the planet's recoil conserves momentum and must depend on 1/m. 

To calculate the trajectory function g{x), we consider nearby particles in the co-orbital zone of 
the planet and more distant particles in the small-angle limit. Planetesimals in the co-orbital zone, 
with \6r\< 2rH, follow almost the same orbit as the plan et but get pushed gently towards and away 



from it on horseshoe orbits ( Dermott &: Murray 



1981 



More distant planetesimals at |(5r|>4rH 
stream by the planet and experience small-angle scattering relative to their Keplerian path. The 
trajectory function in these two cases is 



9{x) 



2x ( I a; I ;S 2; co-orbital), 
-32x~^ (|x|^4; small-angle scattering). 



(3) 



The co-orbital zone result follows from conservation of energy. When the pair's relative speed is 
much greater than the planet's escape velocity at closest a pproach, the small-angle expression for 
larger separations follows from two-body scattering theory (jLin &: Papaloizoulll979l ). 
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To illustrate this scaling property, we consider a set of numerical simulati ons of planetes i- 



mals on circular orbits clo se to a much more massive planet (sim ilar to Fig. 5 inllda et al 



Bromley &: KenvonI (|2006l ) describe our orbit integrator (see also iBromlev &: KenvonI l20ld ) . The 



2000) 



planet and the planetesimal start 180° out of phase on circular orbits at distances a (planet) and 
a + 5r (planetesimal) from the central 1 Mq central star. We measure da as the change in a when 
the planet and the planetesimal complete a single synodic orbit. 

Our results agree very well with the scaling law (Fig. [1]). For three planet masses, scaled 
according to eq. ([2]), the calculated 6a tracks the prediction well for co-orbital particles (|x| < 1.8) 
and distant particles in the small- angle limit (|x| ^4). This scaling law begins to break down when 
the mass of the planet approaches the mass of the central star, but our results lend strong support 
for the "universality" of the trajectory function g{x). 

Orbital separations too distant for co-orbital encounters and too close for small-angle scattering 
encounters lead to chaotic scattering. A formal outer boundary for this limit is 



2\/3rH. 



(4) 



Outside this separation, there is an energy-an gular monientum barrier that prevents chaotic orbit 
crossings for bodies on initially circular orbits (|Gladmanlll993l ). The inner boundary is the edge of 
the co-orbital region; thus, we adopt a chaotic zone with 1.8 rn ^ (5r < 3.5 rn- Fig. [1] shows that 
particles in this region do not follow a simple trajectory function as in eq. 

The trajectory of a particle passing by a planet depends on whether the particle's approach is 
inside or outside the orbit of the planet. This asymmetry is evident in Fig. [H where particles with 
large negative x have smaller \da\ than particles with large positive x. We can correct eq. ^ for 
this property of the orbits using a first order Taylor series expansion: 



da 



9{x) 



(3x^rn dg 
a dx 



m 



(5) 



From numerical simulations, we estimate /3 = 3/8 in the chaotic regime and /? = 9/20 in the small- 
angle limit. Fig. [2] shows 5d in simulations with various planet masses and with planetesimals that 
start inside and outside of the planet's orbit. After normalizing our results using eq. ([5]), these 
traces yield nearly the same "universal" curve g{x). 



2.1. Theoretical migration rates 



To estimate a migration rate from this formalism, we need a relation for the encounter fre- 
quency. Alth ough this frequen cy can vary substantially between consecutive passes of the same 
planetesimal (|Kirsh et al.ll2009l ). a good characteristic number is the inverse of the synodic period. 



1 



3|5r| 



Tsyn 2aT \ Aa J 



1 



55r\ 



(6) 



where T = 27ra^/^/(G'M)^/^ is the orbital period of the planet. The product of this expression and 
eq. ([5]) yields a migration rate per unit mass of planetesimals. 
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Extending this rate to a disk of planetesimals passing by the planet requires a surface density 
distribution for the disk. We adopt a smooth surface density S over an annulus with area 27rr5r. 
Expanding all terms with r = a + 6r in a Taylor series, keeping only first-order terms in 6r, and 
converting to a form with x and dx yields an integral for the migration rate: 



da 
111 



vra^S a 
M T 

1 + 



\x\g{x) dx X 



a dH j]^ ^ 
T, da g dx 4 



rux 



(7) 



where g{x) is from eq. ([3]). The surface density S is often parameterized as a power-law, with 



S(a) 



and n = 1-1.5. We use this form of S throughout, setting n 
unless otherwise specified. 



(8) 



1, oo = 1 AU, and Sq = 30 g cm 



2.1.1. Migration from small-angle scattering 



To understand the implications of eq. ([7]), we consider several simple cases. For distant en- 
counters between a planet and material in a power-law disk at separations > 4, we derive the 
migration rate from eq. ([7]) with g{x) = — and /3 = 9/20: 



da 
It 



327ra2S a 
M T 



/dx 
sgn{x) — 



1 + (2 - n) 



r^ix 



(9) 



where S is evaluated at the planet's position. If a planet lies embedded in a large disk where the 
inner (outer) radii are small (large) compared to the planet's semimajor axis, the first term in 
the square brackets vanishes; the r nigration rate then depends weak ly on planet mass through the 
second term involving rn (see also llda et al.ll2000l : iKirsh et al.ll2009l ). If a planet lies on the inside 
or the outside of the disk, the first term dominates. 

Providing it is no closer than Txing 

from the planet, eq. ([9]) predicts that a disk situated just 
inside or outside a planet's orbit is repulsive. To confirm this behavior numerically, we use an 
A^-body code that evolves massive planets, along with massive planetesi mals that interact with th e 
planets but not with each other (i.e., the disks are not self-gravitating; see lBromlev &: Kenvonll20ld ). 
For example, a 100 M^, n = 1 power-law disk consisting of 2 x 10^ equal-mass particles between 
26.5 AU and 35.5 AU pushes a 0.3 M® planet at 25 AU inward with a speed of 0.012 AU/10 kyr 
(eq. [9]). Simulations of 100 planetary orbits yields 0.0117 ± 0.004 AU/10 kyr. Tests with a disk 
inside the planet's orbit confirm that the planet migrates outward, as expected from eq. ([9]). 

With this formalism, we can consider an idealized example of the migration of a planet nestled 
between two equal-mass annuli of planetesimals. For a surface density S oc a^^ and spacing 
between the planet and each annulus of 6r > rxing, eq. ([9]) predicts a net inward migration. For 
an Earth-mass planet at 25 AU, between two 0.5 AU annuli centered on 23 AU and 27 AU, with 
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50 apiece, the theoretical migration rate is —0.019 AU/10 kyr. Although the planet eventually 
migrates through the gap into the inner disk of planetesimals, our numerical simulation using 
1/600 M0 planetesimals yields a migration rate of —0.025 it 0.002 AU/10 kyr. In this simulation, 
migration leads to more interactions with the inner disk of planetesimals than predicted by the 
analytic theory; still, this numerical result agrees reasonably well with the analytic prediction. 

In these examples, the gap between the planet and the disk spans the chaotic and co-orbital 
zones. In small-angle scatt ering, migration is fairly small, ^ 1-2 AU Myr~^. F or typical growth 
times of ~ 1-3 Myr (e.g., iKenvon &: BromlevI l2006l : iBromlev fc KenvonI I2OIOI ). planets migrate 
through a small fraction of the disk. 



2.1.2. Fast migration 



For a planet embedded in a pl anetesimal disk, the co-orbi tal zone is much more important than 



the small-angle scattering regime (jWard 



1991 



Ida et al.ll2000l ). Over a complete libration period of 



a horseshoe orbit, there is no net migration of a planet responding to a planetesimal. If the planet 
is already moving radially inward (or outward) on a timescale shorter than the libration period, 
the situation is different. The planet can then pull itself along, continually exchanging places with 
the co-orbital material in its path. When this mechanism works, it is efficient and relatively fast. 



Simulations performed with our code and othe r codes (llda et al 



gest that fast migration can be inward or outward. iKirsh et al.l ()2009l ) identify a strong preference 



200nl : lKirsh et al.ll2009f ) sug 



for inward migration. Our calculations confirm this conclusion; more massive planets also seem to 
migrate inward more often than less massive planets. We speculate that inward migration domi- 
nates in most simulations from the gentle inward push of the weakly scattered disk, whose influence 
on a planet increases with m. 

The fast migration rate da{g,st/dt follows from integrating eq. ([7]) over the half of the co-orbital 
zone that a planet traverses. We adopt this half-width as 6r = Xcorn, with Xco = 1.8. However, 
fast migration occurs only if the rate allows a planet to clear the co-orbital zone during the libration 
period of the planetesimal at the zone's edge. Otherwise, the planetesimal orbits back and provides 
a counter-torque before the planet migrates away. Large planetary masses have large co-orbital 
zones that are hard to traverse in a single libration period. Thus, this requirement sets a mass limit 
on fast migration, 

'27ra2EXo 



mfast ~ 4.0 



3/2 



M. 



(10) 



3M 1.8 

In a disk with a surface density of 30 g cm~^ at 1 AU from a solar mass central star, this limit is 
mfast ~ 0.025 Me. 

To estimate the migration rate for planets more massive than mfast, we consider a simple 
model. Fast migration relies on a planet crossing the co-orbital zone, with size 5r ~ rn, within a 
typical synodic period of an orbiting planetesimal. When m > mfast, the co-orbital zone is too large 
for the planet to cross in a single synodic period. Thus, a fraction of the material in the co-orbital 
zone has multiple interactions with the planet, slowing the migration rate. This fraction increases 
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with m; thus, more massive planets migrate more slowly. To quantify this statement, we define 
^cojast as the size of the co-orbital zone for a planet with m = rrifast- Planets with m > mfast 
have larger co-orbital zones, with Xco > Xcojast- For these planets, we assume that planetesimals 
within a distance Xcojast of the planet contribute to migration; co-orbiting planetesimals beyond 
this distance do not contribute. Integrating over this annulus, as in eq. ([7]), and using rjj oc m^/^, 
the attenuation factor scales as m{g,st/n^- For m > mfast, migration rates scale inversely with the 
mass of the plane10. 

Fig. [3] illustrates several nuinerical simulations of fast migration for planets with a broad range 
of masses. Following iKirsh et al.l (|2009l ). each planet lies embedded in a power-law disk extending 
from 14.5 AU to 35.5 AU with E = 1.2 (a/25 AU)~^ g cm~^. We represent the disk with particles 
that are each 1/600*^ of the mass of the planet. To speed up the onset of fast migration, the co- 
orbital zone {6r < rn) is initially clear of particles. We also scale the r.m.s. planetesimal eccentricity 
to keep the same initial e = rn/a for each planet. The nearly identical migration tracks for masses 
below mfast ~ 3 illustrate fast migration at the theoretical rate indicated by the dashed curve. 
When the planet encounters the inner edge of the disk, the rates fall to zero (and sometimes reverse 
sign). More massive planets follow tracks that reflect the 1/m attenuation of the fast migration 
rate for m > mfast- 

Several aspects of protoplanetary disks conspire to set limits on the minimum planet mass for 
fast migration. In a gaseous disk, small planetesimals with St < a are entrained in the gas, where 
St = rpgQ/pCs is the Stokes number, r and p are the radius and mass density of a planetesimal, Pg 



is the local gas density, and Cs is the sound speed (see lYoudin &: Lithwick 1120071 : IChiang &: Youdin 



2010 



Ormel &: Klahri |2010| . and references therein) . Fast migration requires that the Hill radius 



of the planet exceed the scale height, hs, of the planetesimals. Following lYoudin fc Lithwick 
(j2007l ). hs = h min(l, ya/Si), where h is the scale height of the gas and a is the disk viscos- 
ity parameter. Adopting a simple expression for the disk scale height, h = /io(a/l AU)^/^ (e.g.. 



Kenvon &: Hartmann 



1987 



Chiang &: GoldreichI 119971 ) and requiring hg < rn yields a simple ex- 



pression for the minimum mass for fast migration in a gaseous protoplanetary disk: 



k 







0.033 



AU/ ^ 



(11) 



where fst = min(l, (a/^t)^/^). When fgt = 1, eq. ([TT]) yields an ap proximate condition for fast 



(type III) migration t hrough the gaseous disk (e.g., eq. [28l se e also iMasset &: Papaloizou 



2003 



D'Aneelo etaP bood : ICrida et aklbood : IPeplihski et al.l EoOSal . 2008b). When most of the solid 
material is in much larger particles with St ^ 1, lower mass planets undergo fast migration through 
the planetesimals. For 1 km planetesimals with St ~ 10"^ and a = 10~^, fnfast,min ~ 10-6 Me at 
1 AU. 

In a planetesimal disk, particle growth sets another limit on the minimum mass for fast migra- 
tion. During oligarchic growth, leftover planetesimals have typical velocity dispersions, v ~ evesa 



^ Using a different approach, WardI (1991) notes that migration saturates when the planet cannot drift across the 
co-orbital zone in a synodic period (see also Paardekooper &: Papaloizou 2009al ). Our derivation yields the mass 
dependence directly. 
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where Vesc is the escape velocity of the largest oligarch and e ?a S ^/S.^ is the ratio of the sur- 



face ( density in oligarchs to the surface density in planetesimals (e.g., iGoldreich. Lithwick. &: Sari 



2004 : iKenvon &: BromlevI l2008l . 2010). The scale height of the planetesimals is then hg = vVt 



efesc^ ^! leading to a simple expression for the ratio of the scale height to the Hill radius in a disk 
surrounding a solar-type star: 

20pi/6e . (12) 



hs_ 



Thus, planets undergo fast migration through planetesimals only when they contain no more than 
~ 5% (e < 0.05) of the mass in solid material. 



2.1.3. Migration rate summary. 

Here, we summarize the migration rates calculated from eq. ([7]) for fast migration (with the 
reduction factor for large masses), and for a planet embedded in a disk that moves relatively slowly 
through small-angle scattering: 

"dT " ^^'^ ~W Vl~8 / ™i"(l'"^fast/m) - , (13) 

" "s^^"""^ ^ IsmJ JWf' ^^^^ 

where 6R is the physical distance separating the planet and the edge(s) of the disk. The it sign for 
fast migration indicates it can be either inward or outward, at least for small mass planets. In each 
case, we have kept only leading order terms in rn/a and have assumed that the inner and outer 
edges of the disk are far, far away from the planet. 



In a thorough analysis of the orbits of moonlets embedded in Saturn's ring system, ICrida et al 



(|2010l ) derive g{x) in the chaotic regime. Their eq. (39) for the migration rate has the same 
functional form as our eq. (jl5p . including the m}/^ dependence. Because we derive rates in the 
small-scattering l imit, the numerical coefficient in eq. (|15p is a few times smaller than the equivalent 



coefficient in the ICrida et al.l (|2010l ) rate. Considering the differences in the two approaches, the 
agreement in the functional form and the magnitude of the migration rate is encouraging. 

The migration timescale \s t = a/\da/ dt\. For a planet embedded in a power-law disk (eq. [8]), 
the fast migration mode and slow, embedded migration yield 

7"fast(yr) = 2.3 X 10"^ max(l,m/mfast) (16) 
. {5R/af { m 



Wyr) = 3.6x10^^ ) (17) 



5R^ frn ^ "^^^ 



3.6 X lO'' ^— ( — ) (18) 

'"xing 

for fiducial parameters of M = 1 M©, a = 1 AU, and Sq = 30 g cm^^. For other situations, these 
timescales vary as 

r oc a'^'^/^T.-^^M^ (19) 
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where 6=1/2 for fast migration with m < mfast, 6 = 3/2 for embedded migration at fixed 5R/a 
(eq. dZ]), and 6 = 5/6 for embedded migration at fixed dR/vxing- The dependence on these 
parameters is more complicated for attenuated fast migration, with m > m-fast (cf- eq. jlOj). 



Although fast migration is two orders of magnitude faster than the slow mode, embedded 
migration may sometimes dominate. In eqs. (I17p and (llSp . 6R is the distance between the planet 
and the nearby edges of the disk. By construction, 6R > rxing- If the disk is dynamically warm, all 
interactions except for the small-angle scatterings are washed out. In a dynamically col d disk, a 
plane t might make its own gap by scattering away all but the more distant material (e.g., iRafikov 
These two situations are similar to Type I and Type II migration in a gaseous disk (^. 



For many of these migration modes, the power law variation of r with a in eq. I\19\i yields 
an integrable expression for a{t). In all mo des of fast migration arid emb edded migration with a 
constant or slowly varying ratio 5R/a (e.g., [Alexander fc Armitagd l2009l ) . we can adopt da/dt = 
Ca^^ , where 71 = 3/2 — n, and derive a simple expression for the time variation of the semimajor 
axis, 

Ct 



a{t) 



272- 



- + a(0)^/^2 



72 



(20) 



wher e 72 = 2/ (2n — 1 ). Standard models for the minimum mass solar nebula (e.g.. IWeidenschilling 



19771 : iHavashi 1119811 ) have n = 3/2; the semimajor axis of the orbit then contracts linearly with 
time . Radio observation s of young stars are more consistent with n = 1 (e.g.. lAndrews &: Williams 



2007 



Isella et al. II2009I ): a then contracts quadratically with time. Fig. O compares our numerical 



simulations with the analytic results for a{t). 

To apply these results to more realistic disks with planets and planetesimals, we consider several 
examples derived from our planet formation simulations. We begin in §2.2 with disks stirred by 
growing protoplanets, continue in §2.3 with disks declining rapidly in mass, and conclude in §2.4 
with disks containing many growing planets. 



2.2. Stirred disks 



In the planetesimal theory, planets form by accreting smaller objects along their paths. When 
planetesimals are larg e {r > 0.1 km), growing planets gravitationally stir up the orbits of neigh- 



boring planetesimals (|Artvmowica 119971 : iKenvon &: Luul 119981 : iKenvon &: BromlevI |2002| ). Thus, 
growing planets dynamically heat up the disk. In Hill units, with en = ea/rn, ensembles of 
growing planets rarely produce planetesimals with en ~ 5, where e is the orbital eccentricity of 
a planetesimal. During runaway growth, planets grow to masses of roughly 0.001-0.01 M0; the 
r.m.s. eccentricities of planetesimals usually drop from en ~ 100 to en ~ 5. Throughout oligarchic 
growth, e oc £fesc/^X) where v^sc is the escape velocity of the most massive planet, vx is the 
local circular velocity, and e is the ratio of the mass in oligarchs to the mass in plan etesimals 
(e.g., kenvon fc l"^Ii998I : IColdreich. Lithwick. Sarilbooi [Kenvon BromlevI bood l. In Hill 
units, ch oc evesc/rn oc e. Thus, en grows slowly as oligarchs accrete more and more planetesi- 
mals, reaching en ~ 100 during the late stages of oligarchic growth and throughout chaotic growth 



- 10 - 



( Kenvon &i Bromlevll2004l . 2006, 2008, 2010). Inclinations are typically half these valueq 



In a hot disk, interactions be tween a planet and surrounding p lanetesimals weak, slowing the 
migration rate by a factor of (see llda et alJl2000l : iKirsh et alJl2009l ) 



da 
lit 



da 
lit 



eH=0 



[1 + (eH/3)= 



(21) 



To test this prediction, we repeat the calculations for Fig. [3] and vary the r.m.s. value of the initial 
en for planetesimals from unity to 50. For larger values of initial en, migration is undetectable 
in a 10^ yr time frame. For en ~ 10-50, it is a challenge to prohibit particles from the co-orbital 
zone at the start of the simulation. To keep all calculations in this suite on the same footing, we 
allow planetesimals for all initial en to reside in the co-orbital zone. The fraction of particles in 
the co-orbital zone is small; most are not on horseshoe orbits. Still, the onset of the fast migration 
mode is somewhat slow compared to the results in Fig. [3l 

The results of these simulations (Fig. U]) follow the trend expected from eq. (j2ip . Planets 
embedded in a disk of low eccentricity planetesimals migrate rapidly, at a rate that scales inversely 
with the mass of the planet (see also Fig. As we raise the initial en, the migration rate slows. 
For en ~ 50, the migration rate is negligible. As shown by the dashed lines in the upper panel of 
the figure, the reduction in the migration scales roughly as e^^. 



2.3. Eroded disks 



As an individ ual planet migrates through a disk annulus, it disrupts the disk (llda et al.ll2000l : 



Kirsh et al.ll2009l ). In fast migration, a planet tosses material into relatively eccentric orbits in 
random directions. This scattering process reduces the local surface density of planetesimals. 
If the planet moves fast enough to encounter only the unperturbed disk upstream from it, this 
disturbance has little adverse impact on the planet's migration rate. 

For more massive planets, slower migration may lead to a continuous loss of disk material. 
Planets migrate more slowly in less massive disks. To quantify this effect, we let the disk surface 
density within rn of the planet vary as 

S(a,t) = --^S(a,0) (22) 

where t is the time since a planet reaches a narrow annulus of the disk at orbital distance a and e 
describes the efficiency of a planet in scattering disk material. To derive the migration rate for a 
disk with this exponentially decaying surface density, we integrate eq. (j22p over the time it takes 
the planet to migrate a distance rn from eq. ([7]). Thus, we approximate the instantaneous surface 
density by a temporal and spatial average inside an active zone of width O(rH), where e ~ th/o is 
an efficiency factor for clearing the zone of planetesimals. This approximation leads to a non-linear 



•^In a disk where oligarchs grow from collisions with small planetesimals coupled to the gas, e and i are set through 
interactions with the gas instead of stirring by oligarchs (§2.1.2). When the gas produces e > en, migration slows as 
in eq. [21] 
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equation where the average value of S depends on the migration rate. Solving this equation leads 
to a new migration rate: 

2>eri 



da 3erjj 



dt ' 2aT 2ahoT J \ ' ^'^^^ 

where do is the theoretical migration rate of the planet without any time-variation in the disk 
surface density. The sensitivity of this model to the exact form of e is fairly weak. 

When the argument of the logarithm in eq. (j23p is zero, the migration rate vanishes. Thus, 
this expression implies a high mass cut-off, merodo where more massive planets cannot migrate 
through the disk. From eq. ()13p . this limit is 



merode ~ VSvra^Smfast (24) 

/ a \7/4-n/2 

- 0.8 (—) Me. (25) 

If disk erosion is an important process, it begins when a planet reaches roughly an Earth mass at 
1 AU and over 40 at 25 AU. Removal of disk material by scattering is at least enhanced, if 
not entirely enabled, by secondary scattering from nearby planets. Migration in principle, can be 
virtually halted for higher mass planets if they have neighbors that prevent the return of scattered 
planetesimals. 



2.4. Multiple planets in a disk 



In the coagulation paradigm, planets grow hierarchically from smaller planetesimals. As they 
grow, planets almost always have neighbors of comparable mass. During runaway and oligarchic 
growth, a few large objects try to accrete all of the mass in an an nulus. Once these oligarchs contain 
roughly half of the total in ass, they begin to interact chaotically (iGoldreich. Lithwick. &: Sari II2004I : 
Kenyon fc Bromley! l2006l ) . During chaotic growth, planets scatter planetesimals to large en and 
grow by large collisions with other planets (jKenyon &: Bromleyll2006l ). Once chaotic growth begins, 
smooth migration through a sea of planetesimals is impossible. Thus, we consider migration in a 
disk of growing oligarchs which contain less than half of the mass in solid material. 

Planets affect the migration of a neighboring planet in two ways. As a planet migrates through 
a disk, it stirs up the planetesimals along its orbit. After migrating past these excited planetesimals , 
the planet leaves behind a wake of planetesimals with large en (see also lEdgar &: Artymowical2004l : 
Kirsh et al.ll2009l ). This wake is a barrier that prevents other planets from migrating inward from 
larger a. For a planet migrating through a disk with initial en ~ 1, planetesimals left behind have 
typical en ~ 3 — 5. From eq. ()2ip . planets encountering stirred up planetesimals have factor of 2-6 
times smaller migration rates. In addition, planets migrating into a wake require longer periods 
to clear their co-orbital zones of dynamically "hot" planetesimals. As a result of these factors, 
migration ceases. 

Fig. El illustrates this phenomenon. Two planets migrate inward in the fast migration mode; 
the migration of the outer planet stops abruptly when it encounters the wake of planetesimals 
already stirred up by its partner. 
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Migrating planets can also deflect planetesimals that chaotically scatter from a neighboring 
planet. When a planet deflects planetesimals from its neighbor, it prevents the planetesimals from 
returning to the neighbor. The loss of these encounters prevents the neighbor from migrating 
towards the planet. Thus, the two planets recoil from the material that is passed between them 
(see also lMalhotralll993l : Eahn Malhotralll999l ). 



To demonstrate this process, we simulate the migration of a Saturn-mass planet at 10 AU 
embedded in a massive disk (400 between 6 AU and 20 AU, with a power law surface density 



distribution and n = 1 ; see iLevison et al.l 120071 ) . We then vary the mass of a second planet at 



5 AU. As the mass of the inner planet falls from a jovian mass to 30 M®, the sense of migration 
of the outer planet changes from outward (Fig. [6l black curves) to inward (blue curve). Until it 
encounters scattered planetesimals from the outer planet, the inner planet slowly migrates inward 
at the "adjacent" rate from eq. ([9]). Despite its small Hill radius (rn = 0.155 AU), the 30 inner 
planet has a considerable impact on the migration rate of a much more massive outer planet 5 AU 
away. Migration is remarkably fragile. 

To conclude this s ection, we consider migration in a multi-planet syste m. In our simulations of 
planet formation (e.g.. iKenvon &: Bromlevll2006l : iBromlev fc Kenvonll2010l ). planets with masses of 
0.1-1 M^ are often separated by 10-20 mutual Hill radii. To investigate migration in an idealized 
version of these calculations, we simulate the evolution of six 0.5 M^ planets in a disk of plan- 
etesimals extending from 7-35.5 AU. As in the calculations for Fig. [3l the planetesimals have an 
initial surface density distribution S = 1.2 (a/25 AU) 
co-orbital zones are initially fllled with planetesimals. 



^ g cm ^ . Unlike the calculations in Fig. [3l 



Fig. [7] summarizes the main results of these simulations. In a multi-planet system, long-term 
migration rates are small. Initially, each planet clears its co-orbital zone of material in 3 — 6 x 10^ 
yr. Fast migration commences. Eventually, each planet encounters the ensemble of stirred up 
planetesimals left behind by its inward neighbor. Migration stops. In these examples, migration of 
the innermost planet ceases when it reaches the inner edge of the disk. In disks with small inner 
radii, migration of the outer planets ceases well before the inner planet reaches the inner edge of 
the disk. 



2.5. Migration and planet formation 

To place these results in the context of formation scenarios, we consider the growth and mi- 
gration of planets in the planetesimal theory. In this picture, planetesimals ranging in size from ~ 
0.1 km to ~ 100 km condense out of the gaseous disk. Planetesimals collide and merge into larger 
and larger objects. After short periods of orderly and runaway growth, the largest objects enter 
oligarchic growth, where they continue to accrete and to stir up leftover planetesimals. During this 
phase, dynamical friction between planetesimals and oligarchs dominates dynamical interactions 
among oligarchs. Thus, oligarchs remain fairly isolated from one another. Once oligarchs contain 
roughly half of the mass in solid material, their mutual dynamical interactions dominate dynam- 
ical friction with planetesimals. Oligarchy ends. Chaotic gro wth, where oligarchs grow by giant 



impacts and continued accretion of small planetesimals, begins (jGoldreich. Lithwick. k. Sari II2004I : 



- 13 - 



Kenvon &: BromlevI 120061 ) 



During the trans ition from oligarchic to chaotic growth, the 'isolation mass' sets the m ass of 
the largest oligarchs (jLissauen 119871 : iKokubo Idal 119981 : iGoldreich. Lithwick. fc Sari 1 12004 ) . By 
definition, isolated obje cts have smal l dynamical interactio ns; thus, their typical separations are 
~ BvH with B 7-10 (|Lissauejll987l : Ixokubo fc Idalll998l . 2000, 2002). When an object contains 
all of the mass in an annulus of width Brn, it reaches the isolation mass. With mjso = 2TraT,Brti 
and S = Soa-", miso = (27r5Eo)3/2 (3M)-i/2 a3-3"/2. If we adopt a disk with n = 1, Sq = 
10 g cm~^, and B = 7, isolated objects have separations of 2rxing = Brn and lie well outside the 
co-orbital zones of their nearest neighbors. The isolation mass is theri 



mis 



0.07 



10 g cm 



-2 



3/2 



a \3/2 

Tmj) 



1 Mp 







M 



1/2 



(26) 



For each oligarch, the ratio of m-iso to mfast (eq. [ID]) sets the importance, the mode, and 
the timing of migration through a sea of planetesimals. Low mass oligarchs with rrioii < miso 
can migrate, but they cannot migrate freely. The typical radial spacing of low mass oligarchs is 
'i^oii ~ 7(moii/miso)^^^'^H- With r^u <C 7rH, an oligarch migrates only a few rn before it encounters 
the wakes of other oligarchs. Migration then ceases. Once massive oligarchs have m > miso, 
they are free to migrate. However, massive oligarchs also interact chaotically. At 1-10 AU, the 
timescale for chaotic growth is shorter (longer) than the timescale for slow (fast) migration. When 
miso < 'iT'fasti migration is important during chaotic growth. Otherwise, growing oligarchs do not 
migrate through a sea of planetesimals. 

This analysis suggests that migration through an ensemble of planetesimals is rarely important 
within planet-forming disks. From our definitions of miso and mfast (eq- [lO]), the ratio miso/^T-fast = 
3B^/A ~ 37. Although growing oligarchs migrate in the fast mode, they can never migrate very far 
before they encounter the wake of another migrating oligarch. Migration then ceases (Fig. [5] and 
Fig. El). 

There are several plausible exceptions to this conclusion. If collisional damping or gas drag 
reduce the e and i of stirred up planetesimals in the wake of a migrating planet, then another planet 
can migrate through the wake. When the wake consists of large planetesimals with r > 0.1 km, 
however, collisional damping and gas drag are ineffective. Dynamical frictioi i and viscous stirring 
by large planetesimals and small o ligarchs keep particles at large e and i (e.g.. iKenvon &: Luulll998l : 
Goldreich. Lithwick. fc Sari II2004I ). Thus, oligarchs cannot migrate freely through a disk of large 
planetesimals. For smaller p articles, damping inay re duce e and i on timescales comparable to 



the migration timescale (e.g.. iKenyon fc Bromley! l200ll ). Because damping and migration occur on 
similar timescales, closely-spaced oligarchs probably encounter wakes before damping can smooth 
them out. Widely-spaced oligarchs suffer chaotic growth, which keeps planetesimals stirred up 
despite damping. Thus, we conclude that damping does not allow migration in a planetesimal disk. 



Scattering may also lead to effective migration (jLevison et al.l 120071 : [Raymond et al.l l2009al . 



^Our definition is appropriate f or the onset of chaotic growth, when en > 1; when en < 1, the alternative of 
Goldreich. Lithwick. fc Sari I l|2004 ) provides a better measure of the masses of isolated objects. 
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2010)^ 



During chaotic growth, massive planets scatter l ower mass planets farther out in the disk 



2010 



Marzari & Weidenschillins 


2002; 


Veras et al. 


2009; 


Raymond et al. 


2010; 


Bromlev &: Kenvon 



Chatteriee et al.ll2010l ). At large a, oligarchs form slowly (e.g., iKenvon &: BromlevI 120081 . 
2010). Thus, planets formed at small a and scattered to large a may end up in a calm disk 
composed of planetesimals with small ch. Without other oligarchs to impede them, these scattered 
planets can then migrate freely through the outer disk. 

Once chaotic growth ends, any leftover planetesimals can support inward or outward migration. 
For leftovers with large e and i, migration rates are slow. However, outwardly migrating planets 
may reach planetesimals with much lower e and i, enhancing migration rates. Several dynamical 
mod els for the origin of the Solar System rely on migration through a leftover planetesimal disk 



(e.g.. lHahn Sz Malhotra 



1999 



Tsiganis et al.ll2005l ). As planets migrate outward, they may capture 



objects in or bital resonances. This process may yield some dynamical classes of tra ns- Neptunian 



Wvatt et al.ll2005i : .Martin et al. 



objects (e.g., 'Morbid elli et al.l 120081) and dense clumps o f material in debris disks (e.g. IWyatt 



2003 



2007 



Crida et al.ll2009l V 



To conclude this section. Fig. [8] compares the variation of miso and mfast with semimajor axis 
for a plausible disk model. We adopt a disk with S = Ega"^ and Sq = 10 g cm~^ at 1 AU. Here, 
we assume a factor of 3 jump in the surface density of solid material at the snow line, a.^nin 



3 AU . We ignore the likely variation in the position of the snow line with time (jKennedy fc Kenyon 
20081 ). In this disk model, mfast ranges from 0.005 at 1 AU to ~ 1 at 10 AU; miso grows 
from 0.07 at 1 AU to 10 at 10 AU. Based on our simulations, low mass oligarchs with 
m < miso are too closely packed to migrate. Prior to chaotic growth, planetesimals can grow to 
reasonably massive oligarchs. With me^ode > '^iso at all a, disk erosion is also unimportant. Once 
masses reach mjsoi chaotic growth without migration leads to terrestr ial mass planets at 1 AU 



mavmond et al 



2005; 



Kenvon fc Bromley 



at 3-30 AU (jGoldreich. Lithwick. fc Sari II2004I : iBromlev fc Kenvon! l20ld ). 



2006 



Ravmond et al.l l2009bl) and Jupiter mass planets 



3. Relationship to migration in gaseous disks 



To explore whether our results on migration are general, we now consider some analogies be- 
tween migration in gaseous and planetesimal disks. As motivation, numerical simulations demon- 
strate that systems of many oligarchs are likely outcomes of runaway growth in a planetesimal disk. 
Alt hough there are several elegant approaches to the migration of single planets in a gaseous disk 
(see iPapaloizou et al.ll2007l . and references therein) , generalizing these approaches to systems of 20- 
30 (or more) planets with masses comparable to or less than the isolation mass is challenging (see 
Cresswell fc: NelsonI 120061 . 2008). Here, we try to see whether we can apply results for planetesimal 
disks to gaseous disks. 

In the limit of zero v iscosity, equal mass gaseous a nd particle disks provide an identical torque 
on an embedded planet (jGoldreich &: Tremaindll980l ). In both types of disk, local variations in 
density generate the torque. In a particle disk, scattering sets the density structure. In a continuous 
medium, a balance between gravity, pressure, and viscous forces sets the density structure. As the 
viscosity of the medium increases, this structure damps out. In this heuristic picture, planetesimals 
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generate migration efficiently; a very viscous medium cannot generate migration. However, the 
large mass of a gaseous disk gives it an overwhelming advantage over a planetesimal disk. For a 
solar metallicity system, the gaseous disk is roughly 100 times more massive than the disk of solids. 



As predicted by iGoldreich &: Tremaind (|l980l ). our numerical simulations produce coherent 
wakes from orbital resonances close to embedded planets. In planetesimal disks (Fig. [9]), an 
embedded planet scatters planetesimals out of the co-orbital zone into disk regions several rn away 
from the planet. The region of horseshoe orbits is initially filled (as in Fig. [9l left panel); continued 
scattering removes planetesimals from the co-orbital zone (as in Fig. [9l right panel). In both 
panels, a bridge of enhanced planetesimal density connects the high density rings of planetesimals 
lying zb 4-5 rn away from the orbit of the planet. Because the 2:1 resonance lies outside the disk, 
the strongest density enhancements lie at the 3:2, 4:3, and 5:4 resonances. 

To illustrate the time evolution of these structures, the online version of this paper contains 
a movie of a planet migrating from 25 AU to 15 AU in ~ 8 x 10'^ yr. Throughout the movie, the 
planet scatters planetesimals out of its orbit into various resonances. Fig. [10] shows a snapshot 
from the movie. At this point of the evolution, the planet has migrated from 25 AU to ~ 20 AU. At 
25 AU, the original orbit of the planet is nearly devoid of planetesimals. Just outside this orbit, the 
density of planetesimals is somewhat higher than the initial density. Between 20 AU and 25 AU, 
the planet has left behind a sea of stirred up planetesimals, with several density enhancements at 
orbital resonances. At 20 AU, the planet has evacuated planetesimals downstream from its orbit. 
Upstream, planetesimals remain in horseshoe orbits. 

The st ructures in plane t esimal disks are simi l ar to those produced in simulations of gaseous 



disks (e.g.. Bate et al.l 12003; lO'Angelo et al.1 l2003l : Klahr k Kiev! l2006l : IPaardekooper fc Mellema 
2006al ). In all of the simulations of planets within gaseous disks, torques between the planet and 



the disk create local enhancements in the gas density at orbital resonances as well as the bridge of 
material from the planet to the bright rings. For disks with similar surface density distributions 
and planets with similar masses, the derived range of the density enhancements are also similar 
(see, for example IWardlll997l : iD'Angelo &: Lubowll2008l ). Because gaseous disks have some pressure 
support, co-orbital gas lies a small distance, 6rco ~ 0.002 — 0.004a, inside the orbit of a planet (e.g. 
Tanaka et al.ll2002l ). When co-orbital gas lies inside the Hill sphere of the planet (m > 0.03 M^), 
we expect co-orbital gas and planetesimals with similar surface density to exert similar torques on 
a nearby planet. Thus, gaseous and planetesimal disks should produce comparable migration rates 
for planets with m > miso- 



3.1. Migration timescales 

To generalize our migration results to a gaseous disk, we consider the vertical scale height of 
the disk /i as a smoothing length, which sets the minimum size of density features in the disk. In 
Type I migration, this assumption limits the scale and the location of the density wakes that form 
through interaction with a relatively small planet. The largest wakes lie at a distance 6R ^ h ^ r^i. 
Thus, we can use results for small-angle scattering. For gaseous disks, the standard type I rate for 
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an isothermal disk from iTanaka et al.l (|2002l ) is 



dai ,„ „ , , X 27ra^Sgm a 



(27) 



where Tig is the surface density of the gaseous disk at a. Setting 5R ~ h and assuming the 
planetesimals and the gas have the same surface density at the position of the planet, the ratio 
of t he rate from eg. (1141) to this rate is ^/ = 16(2 — n)/3(2.7 + l.ln). Additional features in 
the iTanaka et al.l (j2002l ) derivation, including 3-D effects and corotation resonances, produce the 
different dependence on the surface-density power-law index n. For n = 0.5-1.5, ^/ ~ 2-0.5; 
thus, the rates differ by a small numerical factoi0. Clearly, small-angle-scattering migration in a 
planetesimal disk and Type I migration in a gaseous disk share general properties. 

Type II migration occurs when a relatively massive planet creates a gap in a gaseous disk 
and locks into the disk's viscous flow as a result of a build-up of material at the gap's edges. 
Gravitational torques exerted on the planet by the disk produce inwards migration. Thus, the 
planet responds to the instantaneous density perturbations within the disk. These perturbations 
are strongest at the gap edges, which are several rn away from the planet. The condi tion for gap 



open ing is an elegant inequality between m, h, and the disk viscosity parameter a (jCrida et al 
20061 ): 



3 h 

4 ^ 



+ 50a 



M 



m 



<1 



(28) 

where a = v/h?Q., v is the disk viscosity, and is the angular velocity (e.g., |Pringlelll98lh . This 
constraint has a simple physical interpretation. When the viscosity is small (a — )• 0), the first term 
dominates; planets with Hill radii comparable to the local disk scale height open a gap (see also 
Wardlll997l : iLin &: Papaloizoulll986l . 1993, and references therein). As the viscosity grows (a — )• 1), 
the second term domiiiates; planets with tidal f orces large enough to overcome viscous transport 
open a gap (e.g., IWardlll997l : iBrvden et al.lll999l ). 



From the conditi on for gap opening in t he low viscosity limit and the ITanaka et al.l (|2002l ) 
type I migration rate, iPapaloizou et al.l (|2007l ) derive a simple estimate for the maximum type II 
rate. For planets with m ~ 30-1000 M®, t heir simple estimate ag rees with rates derived from 
detailed numerical simulations (see Fig. 3 of IPapaloizou et al.l 120071 ). Adopting w rn in the 
embedded migration timescale from eq. (jl8p and assuming the same S, we derive a ratio of rates 
'?/7,a=o = 32(2 — n)/3(5.4 + 2.2n) ^ 2.5-0.6 for n = 0.5-1.5. Thus, both approaches yield the same 
scaling with rn and magnitudes consistent to a factor of 2-3. In this limit, irradiation from the 
centr al star sets the scale height of the gas (e.g., iKenvon &: Hartmann 1119871 : IChiang &: Goldreich 
19971 ). Once this scale is set, the mass of the planet establishes the region of the disk that interacts 
most strongly with the planet. In the zero viscosity limit, planetesimals and gas respond to the 
gravity of the planet on the Hill scale, leading to similar timescales. 

In the large viscosity limit, our analogy between planetesimal and gaseous disks breaks down. 
Large viscosity enables a gaseous disk to transport mass inward and angular momentum outward. 



''Using the more recent type I rate from 



D'Angelo fc Lubow 



yields similar results. 
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Viscous transport modifies the universal trajectory function g{x). Identifying g{x) for viscous trans- 
port is beyond the scope of this paper; however, we speculate that substituting proper expressions 
for g(x) and the size of the gap in eq. ([7]) would yield a migration rate reasonably close to published 
type II rates, da/dt ~ a{h/ af'aVl^ where is the angular velocity of the planet at semimajor axis 
a. Successfully applying our approach in the large a limit would link the theories of migration in 
gaseous and planetesimal disks. 

Type III migration is c ompletely analogous to the fast migration mode in planetesimal disks 
(jMasset &: Papaloizoul 120031 . see also eq. [TTj) . Because it regulates how efficiently material is 
transported across a planet's orbital position, viscosity complicates precise comparisons between 
gaseous and particle disks. However, viscosity gener ally makes migr ation in gaseous disks less 
efficient per unit disk mass than in planetesimal disks (jida &: Linll2008l ). As a result, planets with 
masses m uch less than the mass of Saturn are 'safe' from type III migration thro ugh the gaseous 



disk (e.g.. iMasset Papaloizou 



2003 



D'Angelo et al.ll200,'5l : IPephiiski et al.ll2008al . 2008b). 



3.2. Migration with multiple planets 

In a gaseous disk, there are three sources of torque on an embedded planet. Torque from an 
inner spiral density wave and material in the corotation zone produces a net outward migration. 
Material in an outer spiral density wave causes a net inward migration. As the planet migrates 
inward, viscous torques smooth out density perturbations behind the planet. Smoothing occurs on 
a local viscous timescale, which is comparable to the migration rate. 

In a multiple planet system, each planet produces a pair of spiral density waves. Thus, each 
planet feels a torque from the spiral density waves of all planets and the gas in its corotation 
zone. When planets are widely separated, distant spiral density waves contribute little to the 
torque. Widely spaced planets migrate freely. When planets are t ightly packed, many spiral waves 



contribute to the torque. In linear theory (e.g.. lTanaka et al.ll2002l ). multiple torques superpose and 



add to the migration. However, this approach does not address the response of the gaseous disk to 
the time-variable potential of a collection of closely packed planets. The gravitational potential of 
the planets varies on timescales shorter than the viscous timescale, which should wash out spiral 
density waves and reduce migration rates. 

Recent analyses show that the thermod ynamics of the disk is an important factor in setting the 



direc t ion a. nd rate of type I migration (e.g., iPaardekooper k. Papaloizoul l2008l : iPaardekooper et al 



2010 



2011 



In these non-linear calculations, migration of a single planet depends on the vertical 
temperature structure and the relative strength of torques from the corotation zone and the Lind- 
bald resonances. In a multiple planet system, each corotation zone generally lies within a few Hill 
radii of a single planet; thus, closely packed planets may not change the torque from the corota- 
tion zone. Because Lindbald resonances lie many Hill radii away from a planet, they are easily 
perturbed by an ensemble of closely packed planets which change the density and temperature 
structure on timescales shorter than the viscous timescale. Because the disk responds relatively 
slowly to motions of the oligarchs, spiral density waves are probably much weaker in a system 
with many oligarchs than in a system with a few oligarchs. Weaker density waves produce smaller 
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migration rates. By analogy with our simulations of planetesimal disks, we propose that tightly 
packed planets do not migrate. 

To place quantitative constraints on these limits, we compare the locations of the resonances 
that drive migration to the radial spacing of planets. For type I migration, the gaseous disk produces 
the strongest torques at the inner and outer Lindblad resonances, w hich lie at orbital distances 
ScLLR ±2/1/3 from the migrating planet (e.g.. iPapaloizou et al.ll2007l ). For two planets separated 



by r ~ ±4/1/3, their Lindblad resonances overlap. This tight spacing may preclude the elegant 
sprial density wavGS necGSsa^ry for type I migration. Planets separajteci by T^yj^n 

« 2/i have isolated 

Lindblad resonances and can mig rate freely. With h ^ ho (a/1 A\j f/'^ AU (|Chiang k GoldreichI 
1993), this constraint becomes r™„ « 0.06 (a / 1 AVf/'^ AU for ho = 0.03 AU. To convert to Hih 
units, planets with m w have rn ~ O.Ola. Thus, our constraint is 



' mm ^ 6 Tfj I 

m 



1/3 



a \2/7 

1 Au; 



(29) 



Numerical simulations do not yet ad dress constraints on t he ab ility of closely packed planets 
to undergo type I or type III migration. ICresswell &: NelsonI (|2006l . 2008) consider ensembles of 
Earth-mass or larger planets (m > miso) spaced by roughly 5-7 rjj. In their simulations, type I 
migration is briefly interrupted by rapid, chaotic interactions among the planets. Once the planets 
have merged or scattered, type I migration continues. Calculations for systems of lower mass 
planets with m < m-iso do not exist. For now, we assume that ensembles of lower mass planets with 
typical separations smaller than rmin do not migrate. 



3.3. Migration and planet formation 



To establish some constraints on type I migration through a gaseous disk containing an en- 
semble of growing planetesimals, we generalize our discussion of isolated oligarchs from §2.5. As 
planetesimals experience runaway, oligarchic, and chaotic growth, the gaseous disk evolves with 
time. In addition t o viscous evolution, photoev aporation and gas giant planet formation remove 
mass from the disk ([Alexander &: Armitagell2009l'). Observat i ons of young stars suggest t ypical disk 



lifetiines of 1-3 Mvr (IHaisch et al 



2001 



Currie et al. II2009I : Kennedy fc Kenyon II2009I : iMamajek 



2OO9I : I Williams &: Ciezall201ll ). Thus, migration through the gaseous disk ceases after 1-3 Myr. 



The lifetime of the gaseous disk places a roug h lower limit on the masses of planets subject 
to type I migration. For the linear calculations of iTanaka et al.l (j2002l ). the timescale for type I 



migratior i is r oc m ^. T h us, lo wer mass planets migrate more slowly. Setting the timescale in 



eq (9) of Papaloizou et al. (j2007l ) to 3 Myr, we can derive an expression for the masses of planets 
that are safe from type I migration: 



^safe 



0.046 



h/ho 
0.03 

So 



a 



1 AU 

-1 



n+l/14 



1700 g cm 



M 
IM^ 



3/2 



Ma 



(30) 
(31) 
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where we assume a gaseous disk with S = Sq a~". Planets with m < msafc migrate on timescales 
longer than the typical disk lifetime. Fig. [11] compares the variation of msafc and miso with semi- 
major axis. For all a, m\so 3> rrisaic. Low mass oligarchs have long migration timescales; isolated 
objects are not safe from type I migration. 

Despite their lack of safety, many isolated oligarchs are packed too tightly together to migrate. 
To draw this conclusion, we derive the masses of objects with r^m = 7 rn, the typical separations 
of oligarchs at the onset of chaotic growth. Defining rrimin as the mass where r^m = 7 rn 



0.63 



1 AU/ 



\ 6/7 



(32) 



In our picture, oligarchs with m < rrimin are packed too tightly to undergo type I migration. For 
a < 5 AU, rrimin ^ "liso; isolated oligarchs have overlapping Lindblad resonances and are packed 



too closely to migrate. At larger a, 
can migrate. 



< 



3; oligarchs do not have overlapping resonances and 



This result leads to an important conclusion for terrestrial planet formation. If the overlapping 
Lindblad resonances of tightly packed oligarchs at 1 AU do not generate type I migration, chaotic 
growth produces Earth m ass or larger planets on timescales of ~ 10 Myr. In our simulations (e.g., 
Kenyon &: Bromleyl l2006l ) . it takes ~ 3 x 10^ yr (~ 10^ yr) to produce 5-10 (~ 15) oligarchs with 
m ~ 0.002 — 0.004 at 0.85-1.15 AU. These oligarchs are safe from type I migration through the 
gas (Fig. [TT|) . but their low masses allow fast migration through the sea of leftover planetesimals. 
However, growing oligarchs stir planetesimals to en ~ 5. After migrating ~ 0.02-0.03 AU, each 
oligarch encounters planetesimals stirred up by its inner neighbor. Relative to the standard fast 
migration rate, we estimate a factor of 10-100 reduction in the migration rate for each oligarch 
in our 2006 simulations. On the (reduced) migration timescale of > 10^ yr, each oligarch in our 
simulations grows by more than an order of magnitude and begins to interact chaotically with other 
oligarchs. Once chaotic growth begins, oligarchs safely grow into terrestrial planets. 

As chaotic growth ends, several factors probably prevent terrestrial planets from migrating 
through the remnants of the gaseous disk or the sea of leftover planetesimal s . In published sim- 
ulati o ns, leftover planetes i mals have very large e and i ([Ravmond et al.l l2005l : iKenvon Bromley 



2006 



O'Brien et al 



2006 



Ravmond et al.l l2009bl ) : thus, planets can sweep up or scatter the left- 



overs faster than they can migrate through them. For typical disk lifetimes of 1-3 Myr, the reduced 
surface density lowers migration rates through the gas by factors of 10 or more. Migration times 
are then longer than the disk lifetime, saving terrestrial planets from type I migration through the 
gas. 

Formation outcomes for gas giant planets are less clear. In our picture, isolated oligarchs at 
5-10 AU will migrate little through a sea of leftover planetesimals. As chaotic growth begins, these 
objects start to experience type I migration through the gas. The relative importance of chaotic 
growth and migration then depends on several factors. 



1. The masses of leftover planetes imals. Although oligarchs with gaseous atmospheres ac- 
crete small planetesimals rapidly (jlnaba &: Ikoma 1 l2003l : IChambCTsll2006al :l iBromlev &: KenvonI 



2OIOI ). they cannot accrete large planetesimals on timescales shorter than the migration time 
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(e.g., IChambersI l2006bl . 2008). Collisional grinding can reduce the sizes of large planetesi- 
mals, enabling rapid accre tion and the formation of 5-10 cores on very short timescales 
(jKenvon &: BromlevI 12003 ). Thus, rapid core formation depends on the evolution of the size 
distribution of planetesimals during oligarchic and chaotic growth. 

2. The response of the disk to tightly packed oligarchs. When oligarchs are tightly packed, their 
Lindblad resonances overlap. By analogy with our calculations of migration through a sea 
of planetesimals, we speculate that tightly packed oligarchs cannot migrate. However, there 
is no analytic or numerical study of type I migration in gaseous disks with tightly packed 
oligarchs. If tightly packed oligarchs migrate at the 'standard' type I rate, then they migrate 
faster than they grow. If tightly packed oligarchs (and leftover planetesimals) reduce the type 
I migration rate, then they probably grow faster than they migrate. 

3. The response of the disk to planets accreting gas. Once planets reach masse s of 1-10 M, 



they begin to accre t e material froni the disk (jMizuno Ill98d : IStevensonI Il982l : llkoma et al 



200d : iRafikov 1120061 : iHori &: Ikoma Il20ld ). At 5 AU, the nominal migration ti mescale for 10 



planets is ~ 5 x 10^ yr, shorter than the accretion timescale of > 10^ yr ( Pollack et al 



1996 



Bodenheimer et al.ll200d : iKornet et al.ll2002l : iPapaloizou et al.ll2007l : iBromlev Sz Kenvon 



20ld ). However, this planet may not migrate so quickly. When the size of the corotation zone 



is comparable to the disk scale-height, the disk may not be able to launch coherent density 
waves for Type I migration. If corotation torques are important, migration may stall until 
the planet reaches larger masses, forms a gap in the disk, and begins Type II migration 



Masset et al.l (|2006l ). 



4. Discussion 

Fig. [12] summarizes the main conclusions of our analysis. When planets grow in a planetesimal 
disk (left panel), interactions between closely packed oligarchs (m < miso) or between chaotic 
oligarchs (m > niiso) limit migration through a sea of planetesimals. Thus, the building blocks of 
terrestrial planets and ice or gas giant planets are safe from this form of migration. 

In a gaseous disk (right panel), we speculate that low mass planets (m < mmin) are packed too 
closely to undergo type I migration. If this constraint is correct, the building blocks of terrestrial 
planets rarely undergo type I migration. Once they are fully formed, terrestrial planets can migrate 
through the disk. However, the reduced surface density of the disk then limits migration to small 
radial distances. 

Even with these constraints, type I migration is still an issue for the building blocks of gas giant 
planets. At 5-10 AU, gas giant planet formation depends on the relative importance of migration 
and chaotic growth. If chaotic growth dominates, the cores of gas giants can form before they 
migrate. If migration dominates, planets must accrete enough material to begin to accrete gas 
before they migrate into the central star. 

Improving these conclusions requires a better understanding of the transition from oligarchic 
growth to chaotic growth. During the early stages of oligarchic growth, oligarchs are closely packed 
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within a fairly uniform sea of stirred up planetesimals embedded in a fairly uniform gaseous disk. 
As oligarchs grow, they become more and more i solated. As th ey become isolated, oligarchs push 
the excited planetesimals out of their orbits (e.g.. lRafikovll200ll ). This evolution creates two types 
of density perturbations within the disk. 



• At the onset of chaotic growth, closely-packed oligarchs contain roughly 50% of the mass of 
solids. These oligarchs create point-like density enhancements in the surface density distri- 
bution of the solids. 

• Planetesimals contain the other half of the solid material in the disk. Planetesimals tend to 
concentrate in rings between the orbits of the oligarchs. 



Thus, the surface density distribution of the solids is fairly rippled, with planetesimals concentrated 
in the peaks of the ripples and oligarchs orbiting within the troughs of the ripples. 

Current theory addresses the response of the gaseous disk to isolated oligarchs. For a standard 
viscous disk, analytic results and numerical simulations yield reasonably robust solutions to the 
structure of a gaseous disk with an enseinble o f widely spaced oligarchs (e.g.. |Papaloizou et al.ll2007l : 
Cresswell fc: NelsonI 120081 : iLubow &: Idal l20ld ) . Despite many remaining uncertainties in treating 
the (thermo)dynamics of the gas, the eccentricity of the oligarchs, magnetic fields, turbulence, and 
other phenomena, interactions between isolated oligarchs and the disk clearly lead to migration. 

Although the planetesimal theory predicts ensembles of closely packed oligarchs, migration 
theory does not address the structure of the gaseous disk at the onset of chaotic growth. Closely 
packed oligarchs clearly cannot migrate through a planetesimal disk (Fig. [7]). We speculate that 
overlapping Lindblad resonances prevent migration through a gaseous disk. New analytic and 
numerical approaches are required to test this idea. 

Migration theory also does not include the response of the disk to the ensemble of leftover 
planetesimals. Analytic solu tions suggest oligarchs create gaps in the surface density distribution of 
planetesimals (jRafikovll200ll ). Many n umerical simulat i ons show that growing oligarchs push away 



and scatter leftov er planetesimals (e.g. JMalhotra 1119931 : iKokubo Sz Idalll998l : lMorbidelli et al.ll2008 
Kirsh et al.ll2009l . and references therein). With ~ 50% of the solid mass at the onset of chaotic 
growth, structure in the spatial distribution of planetesimals probably leads to density waves within 
the gas. Density waves from individual planetesimals probably have negligible impact on oligarchs 
or planetesimals. However, density waves from the ensemble of planetesimals can interact with 
oligarchs orbiting several th away. It is not clear whether this interaction impacts migration 
significantly; however, including the behavior of planetesimals is necessary for a complete theory 
of migration through a gaseous disk. 

Addressing the response of the disk to closely packed oligarchs and to leftover planetesimals will 
improve our understanding of planet f ormation. Despite our good working knowledge of the growth 
of oligarchs from planetesimals fe.g.. Iwetherill Sz Stewart" 1993 : Kenvon &: Bromley 2008 . 2010), 



the formation of planetesimals (e . g., Youdin 201ol). t 



Goldreich. Lithwick. &: Sari II2004I : iKenvon &: Bromlev 



formed planets within a gaseous disk (e.g., Ilda fc Lin 



le tr ansition from oligarchy to chaos (e.g., 



2006 ) , and the long-term evo l ution of fully- 



20051 : iPapaloizou et al 



2007 



Lubow &: Ida 
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20ld ) are less robust aspects of the theory. Complete numerical simulations of migration with a 



gaseous disk, closely-packed oligarchs, and a sea of leftover planetesimals are beyond the capabilities 
of current computers. Smaller simulations of disks with rings of planetesimals and a few oligarchs 
are possible and would begin to address how planetesimals might change migration rates through 
the disk. 



5. Summary 

We have used analytic results and numerical simulations to explore aspects of migration in 
protostellar disks. 



We derive 'universal' rates for isolated planets migrating rapidly, eq. (|T3|) . or slowly, eq. ([T 
through a disk of planetesimals. When the mass of the planet is much much smaller than 
the mass of the central star, these rates agree with coinprehensive numerical siinulations and 



with rates derived from previou s studies (e.g.. llda et al.ll200Cll : lLevison et al.ll2007l : iKirsh et al 
20091 : iBromley &: Kenyonll201Cl ). We derive an upper limit mfast (eq. (jlOh ) on the mass of a 



rapidly migrating planet. In a disk with surface density S = 30 g cm ^ at a = 1 AU, mfast 



0.025 Mm; for S oc a" 



mfast oc a 



-3 



When m > mfast; fast migration rates are inversely 



proportional to the mass of the planet (Fig. [3|). This result is new. 

• Tests of planets migrating through a disk of stirred up planetesimals verify that r ates scale 
with the eccentricity of background planetesimals in Hill units, e^'^ (Fig- HI see also llda et al 



200d : lKirsh et al.l 120091 '). 



• The strong scaling with ch suggests that planets cannot migrate through the wakes of stirred 
up planetesimals left behind by another migrating planet. Several tests confirm this hypoth- 
esis (Figs. [SHT]). Thus, closely-packed oligarchs do not migrate. This result is also new. 

• When a newly-formed planet migrates or is scattered into a region where planetes imals have 
small ch, this isolated pla net can migrate through a large part of the disk (see also iMalhotra 



19931 : iLevison et al.ll2007l ) 



We use some simple arguments to generalize these results to migration through a gaseous disk. 



Adopting the disk scale height h as the scale for density perturbations in the disk, we show 
that rates for type I, type II (in the zero viscosity limit), and type III migration through 
gaseous disks are similar in magnitude and scaling to rates through planetesimal disks. 

If closely-packed oligarchs migrate as poorly through gaseous disks as they migrate through 
planetesimal disks, we derive limits on the masses of oligarchs that undergo type I migration 
through disks with surface density S = Soa^^. 



Combining these results into a single diagram (Fig. I12p. we conclude that type I migration is an 
important issue during the formation of gas giant planets. The building blocks of these planets are 
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probably safe until they reach the isolation mass (miso; eq. (j26p ). Once their masses exceed miso, 
the migration rate depends on how the gas responds to the mass distribution of smaller oligarchs 
and leftover planetesimals. Addressing this issue requires new analyses. 

For terrestrial planets, we conclude that type I migration is unimportant. Throughout oli- 
garchic and chaotic growth, the building blocks of rocky planets are packed too closely to migrate. 
Once these planets are fully- formed, the surface density of the gas is probably too low to support 
type I migration. Thus, our analysis suggests that standard calculations of terrestrial planet for- 
mation without migration yield robust estimates of the formation timescale and orbital properties 
of terrestrial planets. 

Advice and comments from M. Duncan, M. Geller, D. Kirsh, S. Tremaine, A. Youdin, and an 
anonymous referee greatly improved our presentation. Portions of this project were supported by 
NASA's Astrophysics Theory Program, and the Origin of Solar Systems Program through grant 
NNX10AF35G. 
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Fig. 1. — Derived change in semimajor axis (da) of a planet after an encounter with a planetesimal 
as a function of their initial orbital separation (x) in Hill units. Objects start 180° out of phase on 
circular orbits, with the planet at a = 1 AU from the central star (1 M©); da (scaled by in the 
plot) is the resulting change in a after one synodic period. Planetesimals have masses of 5 x 10~^ 
M®. Colors distinguish planet mass: m = 0.125 (cyan), 1 (black), and 8 (magenta); 
symbol attribute identifies the sense of migration: outward (open) or inward (filled). The scaling 
of da with m agrees with eq. ([2]). The steep curves are theoretical predictions for the co-orbital 
zone (eq. [3]); the horizontal line is from small-angle scattering theory. 
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Fig. 2. — As in Fig. [T] for planetesimal orbits in the chaotic regime. Colors indicate mass for 
0.125 (cyan), 1 (gray/black), and 5 (magenta) planets. The darker (lighter) shades indicate 
interactions with planetesimals that are initially inside (outside) of the planet's orbit. The lower 
panel shows the alignment of these curves after applying the scaling relation in eq. ([5]), which 
transforms da into the "universal" trajectory function, g{x)r}i/m. 
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Fig. 3. — Migration of planets with different masses through a planetesimal disk. The disk is a 
sea o f equal-mass part icles, with S = 1.2(a/25 AU)^^ g cm^^, extending from 14.5 AU to 35.5 AU 
(e.g. iKirsh et al.ll2009l ). Planetesimals have masses of 1/600*'" the mass of the planet and initial 
r.m.s. eccentricity of 1 en- The three lowest mass planets with m < rrifast ~ 3 undergo fast 
migration (heavy dashed curve; eq. [B]), until they "bounce" off the inner edge of the disk. More 
massive planets migrate more slowly (light dashed curves), at a rate that scales with mfast/?7i). 
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Fig. 4. — Planetary migration in disks with S as in Fig. [3] and various values for the initial r.m.s. 
eccentricity. The lower left corner of each panel indicates the mass of the planet. In each panel, the 
initial en varies from 1 to 50 as indicated in the legend of the lower panel (lighter shades correspond 
to smaller initial eccentricity). The dashed curves in the upper panel show predicted rates from 
eq. (EH). 
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Fig. 5. — Disrupted migration of a planet as a result of stirring by its neighbor. The planetesimal 
disk has initial conditions as in Fig. [3l Two 0.5 planets have initial separation of 16 rn- Each 
begins fast migration. The inward motion of the outer planet stops when it encounters the wake 
of excited planetesimals left behind by the inner planet. 



- 33 - 




Fig. 6. — Migration of a Sat urn-mass planet in a planetesimal disk with initial conditions as in 
Fig. 9 of lLevison et al.l (120071 ). The disk extends from 6 AU to 20 AU and has a mass equal to the 
combined mass of Jupiter and Saturn. The Saturn mass planet begins at 10 AU. The inner planet 
starts at 5 AU and has the mass of Jupiter (black) , Saturn (magenta) and 30 M9 (blue) . The cyan 
curve shows the outer planet migrating through the disk in the absence of any other planet. 
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Fig. 7. — Migration of multiple planets in a planetesimal disk with initial conditions as in Fig. [3] 
except that the inner edge of the disk is at 5 AU. Six 0.5 planets separated by 16 rn migrate 
inward until they encounter the wakes of their inner neighbor. Once the innermost planet reaches 
~ 10 AU, its mass exceeds mfast- Because mfast oc a^/^, its migration rate then slows dramatically. 
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Fig. 8. — Growth and migration modes in a planetesimal disk. The heavy sohd curve and the 
shaded region indicate the variation of the isolation mass {rriiso) with semimajor axis; planets with 
m < miso {m > rriiso) undergo oHgarchic (chaotic) growth. The dot-dashed curve indicates the 
variation of rrifast with a; planets with m < mfast (fn > mfast) undergo fast (slow) migration. Until 
planets reach miso, they are tightly packed and unable to migrate large distances through the disk. 
Once they have m > rriiso, they are free to migrate in the slow mode. As planets grow larger than 
miso, their migration may be slowed by disk erosion, as indicated by the dotted line (rrierode)- 
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Fig. 9. — Density wakes in planetesimal disks with an embedded planet. Both images are in a 
frame rotating with the planet, which has a mass of 16 and a semimajor axis of 25 AU. In 
the left panel, the corotation zone contains planetesimals; in the right panel, the corotation zone 
is empty. From the inner edge of the disk at ~ 15 AU to the outer edge at ~ 35 AU, the images 
show the local planetesimal density - averaged over 1 kyr - relative to the initial surface density, 
S(a) = 30 g cm~^ (a/1 AU). In the lower right corner of each image, the scale shows the linear 
map of density to color. The full range of the color map is a factor of two in the local mean density. 



Fig. 10. — As in Fig. [9] for a planet experiencing fast migration. This image is a snapshot from a 
simulation of a 1 planet, available in the on-line version of the Journal. The simulation shows 
density structures after the planet has moved several AU inward from its initial orbital distance 
at 25 AU. When the planet is in fast migration mode, the upstream corotation zone is filled; the 
downstream region is relatively empty. 



- 38 - 




1 10 
a (AU) 

Fig. 11. — As in Fig. [8] for a gaseous disk. Curves for isolation mass are also from Fig. [51 Planets 
with m < uisafe (dot-dashed line) migrate on timescales longer than the lifetime of the gaseous 
disk. Before they migrate significantly, the gas disperses. Planets with m < rrimin (dashed line) are 
packed too closely to migrate through the gaseous disk. Terrestrial planets likely undergo chaotic 
growth before they are able to migrate. The cores of gas giant planets start to migrate as they 
begin chaotic growth. 
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Fig. 12. — Migration in gaseous and planetesimal disks. In a planetesimal disk (left panel), planets 
with m < m-iso are packed too closely to migrate. When m > m[so, chaotic growth dominates 
migration. In a gaseous disk (right panel), planets are spaced too closely to migrate when m < 
iTT'min- Once m > rrimin, planets grow chaotically as they migrate. The relative importance of 
chaotic growth and migration probably depends on the response of the disk to smaller oligarchs 
and leftover planetesimals. 



